First Principle Simulations of Heavy Fermion Cerium Compounds Based on the 

Kondo Lattice 
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the hybridization function between the localized 4f or- 
bitals and the s,p,d-conduction bands calculated by the 
first-principle electronic structure calculation based on 
the local-density approximation (LDA) with Hubbard 
I type of the self-energy for the f electrons. Then 
the Kondo coupling is defined via the Schrieffer- Wolff 
transformation 14| , and the KLM is solved with the new 
efficient CT-QMC Kondo impurity solver [15( combined 
with DMFT. 

Now we define the realistic KLM Hamiltonian for a 
given Cerium compound. The general Coqblin— Schrieffer 
Hamiltonian is the following. 
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We propose a new framework for first-principle calculations of heavy-fermion materials. These 
are described in terms of the Kondo lattice Hamiltonian with the parameters extracted from a 
realistic density functional based calculation which is then solved using continuous-time quantum 
Monte Carlo method and dynamical mean field theory. As an example, we show our results for 
the Neel temperatures of Cerium-122 compounds (CeX 2 Si2 with X=Ru, Rh, Pd, Cu, Ag, and Au) 
where the general trend around the magnetic quantum critical point is successfully reproduced. Our 
results are organized on a universal Doniach phase diagram in a semi-quantitative way. 
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First-principle description of heavy-fermion materials 
has been a challenging problem for a long time. The diffi- 
culty arises from the dual nature of the electrons between 
localization and itinerancy due to the large Coulomb re- 
pulsion energy U at each site of the lattice. Here the rel- 
atively well-localized f-electrons interact with the itiner- 
ant s,p,d-electrons that form the conduction band. The 
heavy-fermion systems are generally described as the An- 
derson impurity model in the dilute limit [l[ or the An- 
derson lattice model (ALM) in the dense limit, and the 
first-principle description ofthem has been done by sev- 
eral authors [2j, y, U, [5[ • With the development of dynami- 
cal mean field theory (DMFT) Q and novel continuous- 
time quantum Monte Carlo (CT-QMC) solvers for the 
impurity problem @, B , the ALM description has be- 
come quite successful except for a very low temperature 
range. The numerically exact treatment of the Ander- 
son impurity problem is still very expensive if the tem- 
perature range of 0(1) K is to be reached. Thus, the 
first-principle description of strongly-correlated materi- 
als around the quantum critical point (QCP) [lfj, which 
has recently been attracting a lot of research interest is 
yet to be solved. 

Here we attack the problem using the Kondo lattice 
model (KLM)trying to focus on the low-energy physics 
of the ALM. We show how this new approach works for a 
archetypical family of so called Cerium 122 compounds, 
CeX 2 Si 2 (X=Ru, Rh, Pd, Cu, Ag, and Au), which 
has been one of the most extensively studied strongly- 
correlated materials since the discovery of heavy-fermion 
superconductor CeCu2Si2 [H|]. Strictly speaking, we deal 
with the Coqblin-Schrieffer model [12J with full 14 fold 
degenerate f-shell but effectively the degeneracy is low- 
ered due to the spin-orbit and crystal-field splittings. 
With the localized Kondo-impurity picture we can save 
the amount of the degrees of freedom in our model by 
eliminating charge fluctuations, and we can reach much 
lower temperature range as compared to the ALM sim- 
ulations. The conduction band in the model is given by 
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Here e k is the conduction band, Jk is the Kondo coupling, 
^splitting 1S t ne crystal and spin-orbital field, cu and fi a 
are the annihilation operators for the conduction and 4f 
electrons, respectively, with the orbital a on the lattice 
site i. To solve this Hamiltonian we first need to define 
Jk and the conduction electron Green function. For this 
we perform the first principle DFT calculation within the 
local density approximation for s,p,d electrons plus the 
Hubbard I approximation for the f electrons based on 
the full- pot ential linearized muffin-tin orbitals (LMTOJ 
method [l3l | and calculate the hybridization function [3] 
3A Q (e) = nJ2 k \V ak \ 2 S(e-e k ) ^ n\V a \ 2 p(e) where V ak is 
the hybridization matrix element and p(e) is the density 
of states of the conduction electrons at energy e which 
we measure from the Fermi energy. We use experimental 
lattice parameters for all materials that we study. 

The calculated Tr$A(e)/(TrN F ) = 

[VWllEl'i 84 "^) is shown in Fi g- HI for sev- 
eral representative CeX2Si2 materials with X=Ru, Rh, 
Pd, and Ag. Here Np = 14 is the total number of 
degeneracy and the trace of 3A is taken over all of 
Np states. We note that SA shows strong frequency 
dependence therefore in order to define Jk an averaging 
over some frequency intervals needs to be performed. 
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FIG. 1: (Color online) The hybridization function 
Tr3A(e)/14 between the conduction band and the 4/- 
electrons calculated by LDA + Hubbard I for CeX2Si2 with 
X=Ru, Rh, Pd, and Ag. The origin of the energy is set to be 
the Fermi level. 



The Kondo coupling Jk is denned by the Schrieffer- 
Wolff transformation [14j, [l6| as follows 
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Here ef is the location of the energy level of 4f orbital, and 
U c tf = U — Jnund is the effective on-site Coulomb repul- 
sion taking into account an effective Hund coupling JHund 
that works in the virtual f 2 state. We set ef — — 2.5[eV] 
and U = 5[eV] which is a typical value for Cerium com- 
pounds. The Hund coupling Jnund is explored around a 
realistic value 1 eV as is explained below. In the present 
formulation, U c g incorporates all of the possible multiplet 
effects in the virtual f 2 states and some systematic error 
comes in from the setting of this value, but it is small 
enough to see the general trend between the materials in 
the realistic Doniach phase diagram that is obtained in 
Fig. 2] in the end. Here we have a band cutoff -D cu toff set 
to be 5 [eV] which is large enough to make a universal 
description of the low-energy physics 17 1. 

The portion of the conduction electron Green func- 
tion G a (e) which has non-zero hybridization with the 
f-electrons is also proportional to A Q (e). We define the 
normalized and Hilbert-transformed G a {iui) as follows 
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The Eqs. J2]), ([3]) provide necessary inputs which are 
plugged into the CT-QMC and solved with DMFT self- 
consistency loop. The details of the CT-QMC algorithm 
for the Coqblin-Schrieffer model are given in [IH . These 
definitions for the realistic model are designed in such 
a way that it becomes exact in the limit of constant hy- 



bridization with the relevant quantity Np Jkp(O) that de- 
termines the behavior of the KLM [15| . 

The LDA results for TVp Jxp(O) for the target materials 
are given in Table |TJ The level splittings A a appeared in 
[T]are implemented as the difference of the positions of ef 's 
which are used in the update probability as is described 
in Ref. [lij]. These level splittings are taken from the 
literature and summarized in Table fl] We checked that 
our results for the Neel temperatures are robust against 
small changes of factor of O(l) on the level splittings. 
These A Q reduce the effective degeneracy close to Nf — 
2 [3]. Thus we call our model "realistic Kondo" lattice 
instead of the Coqblin-Schrieffer lattice even though we 
are actually doing the multi-orbital model. 

We apply the above framework for the KLM descrip- 
tion of CeX 2 Si 2 with X=Ru, Rh, Pd, Cu, Au, Ag. 
We do the following analyses with several settings of 

U c s — U — Jnund f° r < Jnund ~ 1 eV for each of 
the material. For a given material and given parameter 
set, we determine the Neel temperature by looking at the 
temperature dependence of staggered susceptibility and 
locating at which temperature it diverges. Here we fol- 
low the formalism of DMFT for the localized f-electron 
systems as given in [l9| and use the same method as 
was utilized for model calculations in [2(|. Regarding 
the realistic input of the Green's function as is depicted 
in Fig. [1] , we make an approximation in the calcula- 
tion of the staggered magnetic susceptibility for the 4f- 
electrons by decoupling the two-particle density of states 
p( e ii e 2) = S(ei +£2)p( e i) a s if there is a nesting property 
which becomes exact when the 4f-electrons are on a hy- 
percubic lattice. Thus the tendency to the antiferromag- 
netic order would be overestimated in addition to having 
the infinite-dimensional nature in the DMFT solution to 
the lattice problem. The data specific to CeRh2Si2 with 
which we determine the Neel temperatures for several set- 
tings of JHund are shown in Fig. [5] In this way for each 
of the material we look at the magnetic phase transitions 
for several Jk's by varying corresponding JHund 's. 

As was first discussed by Doniach 2l|, |22j and subse- 
quently by many authors, Kondo lattices have two repre- 
sentative energy scales, namely the magnetic ordering en- 
ergy that is proportional to ( Jk/-/Vf) 2 p(0) and the Kondo 
screening energy which behaves like exp (— 1/N-p Jkp(O)). 
For small Jk's the former wins but as Jk becomes larger 
the exponential growth of the latter dominates at some 
point. Thus a given system can realize in either magnet- 
ically ordered phase or non-magnetic Kondo-screened 
phase. Between these two phases at zero temperature 
there is thought to be a QCP. We explore this Doniach 
phase diagram for each material and find the material- 
specific QCP. We take the data with the setting Jnund = 
0.94 eV as our realistic result for each material as this 
strength of the Hund coupling is close to the realistic 
value and also gives reasonable trend over all materials 
in the family. Thus obtained Doniach phase diagrams for 
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TABLE I: Inputs (given by LDA and experiments in the literature) and outputs for each material. A sp i n -orbit is spin-orbit 
splitting between j — 5/2 and j = 7/2 states. Crystal field splitting of j — 5/2 state produces three doublets with energy 
E <E 1 < E 2 with E 1 -E = Al plitting and E 2 - E = A s 2 plitting . J Hu nd ~ 1 eV is used. The Aj plit «^, A s 2 plitting , AJ^, ut 
are given in meV unit and Tk in K. Similar values for the level splittings were used in a recent work [371 ]. 
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FIG. 2: (Color online) Determination of Neel temperatures 
for several settings of the Hund coupling JHund- We see that 
the Neel temperature disappears at the point where JHund 
takes the value between 0.96 and 0.98 eV. 



all of the materials are shown in Fig. [31 

Now we can plot all of the six materials CeX2Si2 
(X=Ru,Rh,Pd,Cu,Ag,Au) on the universal Doniach 
phase diagram in the same spirit as was done by Endstra 
et al. in 1993 [36| but with a different horizontal axis. In 
the material-specific Doniach phase diagrams in Fig. [3J 
we see that the locations of QCP's on the Np Jkp(0) are 
not actually universal (23[. So we measure the distance 
between the QCP and the material's realistic location 
on the horizontal axis, NfJkp{0) — NfJkp(0)\qcp, an d 
plot the Neel temperatures with respect to the value 
of the distance to QCP defined as t = [NfJkp(0) — 
-/Vf Jkp(0)|qcp]/-^f Jkp(0)|qcp- The result is shown in 
Fig. 2J We expect a systematic error bar in the estima- 
tion of the value on the horizontal axis especially around 
the QCP but these possible systematic errors are small 
enough to discern the locations of CeX 2 Si2 (X=Ag,Au) 
and CeRli2Si2. The antiferromagnet CeRli2Si2 and the 
paramagnet CeCu2Si2 are mixed up within the present 
level of resolution which is manifested by the result that 
finite Neel temperature is plotted for CeCu2Si 2 . How- 
ever, this is actually caused by the proximity of this ma- 
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FIG. 3: (Color online) Material specific Doniach phase dia- 
grams for CeX 2 Si 2 with (a) X=Ru,Cu, and Rh and with (b) 
X=Pd, Au, and Ag. The realistic results with JHund = 0.94 
eV are marked with the larger plot symbols. 



terial to its QCP as can be seen from Fig. [3J . So our nu- 
merical result is consistent with the experimental result 
that CeCu2Si2 is a superconductor at ambient pressure 
and is thought to be close to the QCP. 

We note that the valence fluctuations which we ignored 
in our simulation could be important in the realization of 
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FIG. 4: (Color online) Universal Doniach phase diagram for 
the material family of CeX2 Si2. The horizontal axis is defined 
as follows: t = [N F J K p(0) - N F J K p(0)| QC p]/iVp J K p(0)| Q cp. 
The line is guide for the eye. We plot CeCu2Si2 with an aster- 
isk mark to note as it apparently has a finite Neel temperature 
but that just reflects the result that CeCu2Si2 is very close to 
QCP as can be seen in Fig. [3] 



the non-magnetic ground state. Indeed it is known that 
there are some valence fluctuations in CeCu2Si2 24 1 and 
CeRu2Si2 jig ]. This might make the possible systematic 
error relatively larger on the right-hand side of our phase 



diagram [26|, |27J. Nevertheless at the present level of 
description we believe that the realistic KLM works be- 
cause the number of 4f electrons in Cerium ion is still very 
close to one 



24, 



25]. At least for the impurity problem 
the convergence to the Kondo impurity picture in large 
kf |/(jp(0)y 2 ) limit of the Anderson model was discussed 
exactly |28[. Careful comparison between the Anderson 
lattice and the Kondo lattice regarding the valence fluc- 
tuation issue is interesting, especially for CeCu2Si2, and 
further work is ongoing in this direction. 
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